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In the independent-particle model, the nuclear level density is determined from the neutron and 
proton single-particle level densities. The single-particle level density for the positive-energy con- 
tinuum levels is important at high excitation energies for stable nuclei and at all excitation energies 
for nuclei near the drip lines. This single-particle level density is subdivided into compound-nucleus 
and gas components. Two methods were considered for this subdivision. First in the subtraction 
method, the single-particle level density is determined from the scattering phase shifts. In the 
Gamov method, only the narrow Gamov states or resonances are included. The level densities cal- 
culated with these two methods are similar, both can be approximated by the backshifted Fermi-gas 
expression with level-density parameters that are dependent on A, but with very little dependence 
on the neutron or proton richness of the nucleus. However, a small decrease in the level-density 
parameter was predicted for some nuclei very close to the drip lines. The largest difference between 
the calculations using the two methods was the deformation dependence on the level density. The 
Gamov method predicts a very strong peaking of the level density at sphericity for high excitation 
energies. This leads to a suppression of deformed configurations and, consequently, the fission rate 
predicted by the statistical model is reduced in the Gamov method. 
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I. INTRODUCTION 

The nuclear level density p is an essential ingredient in 
calculating the statistical decay of a compound nucleus 
(CN) by particle evaporation, gamma-ray emission, or 
fission. The statistical model has widespread use in nu- 
clear physics and applied research. All told in these ar- 
eas, knowledge of the level density is needed at low and 
high excitation energies, with small and large compound- 
nucleus spins, and for the full range of Z and N from 
around the f3 line of stability out towards the drip lines. 
For example, the cross sections for neutron capture on 
nuclei close to the neutron and proton drip lines are of 
interested in r and rp nucleosynthesis calculations. If 
the excitation energy in these reactions is sufficient, the 
statistical model is used to determined the n and 7 de- 
cay rates of the fused system. In such applications, level 
densities are required for nuclei with extreme n-p asym- 
meties. For many of these nuclei, it will not be possible 
to measure the level density even with proposed radioac- 
tive beam facilities. Clearly a good understanding of the 
dependence of the level density on the n-p asymmetry is 
required to extrapolate to these systems. Even for less 
exotic compound nuclei closer to the f3 line of stability, an 
asymmetry dependence can have important consequences 
on the n-p asymmetry of the evaporation residues[J. For 
fission decay, the deformation dependence of the level 
density is also needed. Therefore it is important to know 
the excitation-energy, asymmetry, and deformation de- 
pendencies of the level density over most regions of the 
chart of nuclides. 

A complete understanding of the nuclear level density 
requires consideration of the many-body nature of the nu- 
cleus. However, the independent-particle model provides 



a useful reference to start with. It also permits a rapid 
survey of level-density dependencies over many regions of 
the chart of nuclides and gives insight into how different 
nuclear-structure effects modify the level density. Many- 
body effects such as the effective-nucleon masses and col- 
lective enhancement due to rotational and vibrational 
collective modes can be incorporated in a phenomenologi- 
cal way onto the independent-particle model0,|j| . In this 
paper, the excitation-energy, n-p asymmetry, and defor- 
mation dependencies of the level density are investigated 
within the framework of the independent-particle model. 
Specifically, the role of the continuum of positive-energy 
single-particle states is studied. For systems around the 
f3 line of stability these states are populated significantly 
only at large excitation energies. However for systems 
closer to the drip lines, where either the neutron or pro- 
ton separation energy is small, these states can be popu- 
lated significantly even at low excitation energies. It has 
been suggested that the contributions from these contin- 
uum states may lead to a n-p asymmetry dependence of 
level density 4] . The manner in which these states influ- 
ence the deformation and excitation-energy dependencies 
of the level density will also be investigated. 

Before further discussion of the level density, it is useful 
to first consider the largest excitation energies for which 
it is meaningful to apply the statistical model. The CN 
is a system of nucleons which is equilibrated in its single- 
particle degrees of freedom and thus has a long lifetime 
compared to the timescale of single-particle motion. As 
such, compound-nucleus decay is a rare process, i.e., the 
typically energy fluctuation of a nucleon, which is of order 
of the temperature T, does not lead to the emission of 
that particle. Thus, the regime of applicability is T < 
ETost where Ef™ t is the minimum of E 7 c l ost and E% oat , 
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the energetic costs of emitting a neutron or a proton, 
respectively. For neutrons, the cost is just the neutron 
separation energy E£ ost =E™ , while for protons the cost 
also includes the Coulomb barrier E^ ost =EP ep + V cou i- 
The decay width for protons or neutrons is roughly^ 



T = exp (-Ecoat/t) 

7re 



(1) 



where e —h 2 /2mR 2 , 1/t — dp/dE* is the nuclear tem- 
perature (t « T), m is the nucleon mass, and R is the 
nuclear radius. As t approaches E^f t , the decay width 
for either proton or neutron evaporation becomes quite 
large. For the statistical model to be applicable, the to- 
tal decay width must be small compared to the spread- 
ing width which determines the time scale for the CN to 
equilibrate. 

The order of this paper is as follows. A brief review 
of the level density in the independent-particle model is 
given in Sec. [H] Subsequently, two methods to include 
the positive-energy states are considered in Sec. lllll Fol- 
lowing this in Sec. lIVI the details of the coupled-channels 
calculation of the single-particle level densities are given. 
The determination of the deformation energy is discussed 
in Sec. I VII and calculated level densities are presented in 
Sec. IVIII Finally in Sees. IVIIII and IIXI a discussion of 
the results and the conclusions of this work are made. 



II. LEVEL DENSITY IN THE 
INDEPENDENT-PARTICLE MODEL 

The independent-particle model starts with sets of 
single-particle levels for both neutrons and protons. The 
determination of the nuclear level density is essentially 
a combinatorial problem, i.e., how many ways can these 
single-particle levels be occupied to give the desired total 
excitation energy. However, the enumeration of all the 
single-particle configurations can be avoided. Instead, 
the Laplace transform Z(a, (3) of the level density is more 
easily calculated when the Lagrange multipliers a and (3 
are introduced to constraint the total number of parti- 
cles and the total energy. The inverse transform can be 
obtained from the saddle-point approximation to give a 
formula for the level density which is continuous in exci- 
tation energy E* . For simplicity at this point, consider 
only one type of nucleon with single-particle levels £j, the 
level density is then@ 



p(E*) 



exp S 



(2) 



where S = (3E - aA + \nZ(a,(3). The values of 
the Lagrange multipliers are determined by the saddle- 
point condition dS/d(3=dS/da=0. Now the average 
occupancy of a single-particle level is given by Si = 
1/ [1 + exp (/fej — a)]. Thus the saddle-point condition 
can be expressed in terms of the conservation of nucleon 



number A and total energy E = E gs + E* (ground-state 
+ excitation energy) by 



da 



E — E gs + E* 



d\nZ 

8(3 



$>i/i- 



(3) 
(4) 



At the saddle point, the quantities Z, D, and S are now 
In Z = ^ In [1 + exp (a - /3e,)] , (5) 

i 

d 2 In Z d 2 lnZ 



D 



do 2 dad{3 
5 2 In Z d 2 lnZ 



and 



Here 



dpda d/3 2 



-/( In Si- (1-/0 In (1-/0 



(6) 



(7) 



(8) 



Although this formula can be derived without recourse 
to statistical mechanics, Bethe 7] realized there is a close 
analogy to the problem of a Fermi gas in contact with a 
heat bath of temperature T—1//3 and with chemical po- 
tential p=a//3. In this analogy, Z is the grand partition 
function, S is the entropy and thus 1/T = dS/dE*. 

If the single-particle level-density g(e) = J2i H £ ~ e 
is constant (at least in the vicinity of e=p), then Eq. [21 
can be reduced to the well known Bethe or Fermi-gas 
expression [T 



p(E*) 
S 



exp S 



V48S* 

2\foE* = ciT 2 , 



(9) 
(10) 



where a—^-g(p,) is the level-density parameter. For 
a two-component Fermi gas, the level density pa- 
rameter will have contributions from each component 

a>=\- [pnO-'n) + 9p{^p)]- Experimentally, level-density 
parameters exhibit strong shell corrections at low exci- 
tation energies. However apart from this, the average 
value of the level-density parameter is often assumed to 
depend only linearly on A with no dependence on the n-p 
asymmetry. 

To gauge the temperatures for which this formula 
should be applied, the functions / (Fermi function) and 
s, which are needed to determined the total energy and 
entropy (Eqs. 01 and |7J|, are plotted in Fig. Q] verses 
f3e — a={e — p) /T. The Fermi function /, giving the 
average level occupancy, changes in value from 90% to 
10% over an interval Ae=4.4T centered around p. The 
function s is Gaussian-like with a full width half maxi- 
mum (FWHM) of 4.2T, however the tails of the function 
falls off much slower than a Gaussian. The Fermi-gas 
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FIG. 1: (Color online) The evolution with (e - p)/T of the 
functions / and s (Eq. used to determine the total energy 
and entropy. The vertical dashed line indicates the location 
of the e=0, when T=E sep /2. 



formula thus assumes the single-particle level density g 
is constant at least over an interval ±2T around p. How- 
ever because s falls off so slowly, the contribution to the 
entropy from levels at smaller and larger energies are not 
insignificant. Therefore at large temperatures, how use- 
ful is the Fermi-gas formula when g is not constant? At 
low temperatures by expanding Eqs.|321 and [3 as func- 
tions of T, the entropy with its lowest-order correction 
becomes S — V a'E* where 



5.9(Mo) 3 



E* 



(11) 



fiQ is the chemical potential at T=0, and g' and g" are 
the first and second derivatives of g. Thus the level- 
density parameter can be replaced by an effective level- 
density parameter a' which is excitation-energy depen- 
dent. Higher-order corrections will be needed at larger 
temperatures. 

What about the role of positive energy states? Con- 
sider a system where T =E sep /2. In Fig. ^for this case, 
e=0 would correspond to the vertical dashed line and 
clearly, because s decreases so slowly with e, the positive- 
energy states (those beyond the dash line) are important 
for the entropy. 



III. CONTINUUM SINGLE-PARTICLE LEVEL 
DENSITY 

In order to understand the role of positive-energy 
states in calculating the level density, we need a pre- 
scription for deciding which of these positive-energy lev- 
els belong to the CN. If there is no bounding volume 
enclosing the nucleus, the single-particle level density of 
these states is infinite. However not all of these states are 
considered to be associated with the CN. To better un- 
derstand the selection of positive-energy states belonging 



to the CN, one can consider the problem of nucleon evap- 
oration. This is often dealt with by placing the CN in a 
box whose volume is large compared to the nuclear vol- 
ume. Call this state i. The box volume V can eventually 
be expanded and allowed to approach infinity. If the CN 
decays by the emission of a nucleon with kinetic energy 
£ to state /, then as we have a bounding box, the nucleon 
cannot escape and will eventually be reabsorbed by the 
daughter nucleus leading us back to state i. Therefore by 
the general principle of detailed balance, the transition 
probability Pi / from state % to / is related to that of the 
inverse process by 



Pi p if = Pf p fi 



(12) 



where pi and pj are the density of states for i and /, 
respectively. For state i, the density of states is just the 
level density of the compound nucleus pi — pcn(E*). 
While for state /, both the level density of the daughter 
nucleus pa and the phase space g ev of the evaporated 
particle contribute, .i.e., 



p f = p d {E* - E sep - e)g ev (e) de. 



(13) 



Now 



9ev(e) 



= (2s + 1) 



(2m) 3/2 V^e~ 
4tt 2 ?i 3 



(14) 



where p, m and s are the evaporated nucleon's momen- 
tum, mass and spin, respectively. Here g ev is deter- 
mined for an empty box in the semiclassical limit. This 
should be appropriate as the box volume is large. Be- 
cause Pfi — v (Ji nv (e) /V, then 



Pa = 



T (e) de 

h 

(2s + l)m 



-EOi 



w (e) 



p d {E* - E s 



PCN (E*) 



tfe(15) 



where v is the nucleon velocity and Ui nv is the inverse or 
absorption cross section. 

In this derivation of the Weisskopf evaporation for- 
mula, the single-particle level density contributes to both 
Pcn and the phase space of the evaporated nucleon g ev . 
Thus for a given nuclear mean-field potential surrounded 
by a bounding box, the total single-particle level density 
will be subdivided; g to t (e) = gcN (e) + g gas (e) where 
gcN is the single-particle level density used to calculate 
the compound-nucleus level density and the remaining 
level density g gas is associated with a gas of evaporated 
particles. Thus g gas ~ gev, where g ev is the single- 
particle level density for the empty bounding box, i.e., 
without the nuclear mean- field potential (Eq. I14f> . As 
the box volume is chosen to be much larger than the 
nuclear volume, then we also find gtot ~ g e v, though of 
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course gtot ^ 9 gas- With such a subdivision of gtot, 
nucleoli number in the box can be subdivided, i.e., 



the 



Atot = I gtot (e) de = gcN (e) de + g gas (e) cfe 



.4 



CjV 



A r 



(16) 



Similarly B 4ot = Ecn + E gas and 5tot = Scn + S gas - 
For a given temperature, the chemical potential \x is 
constrained so that Acn is the constant value appro- 
priate for the CN. Thus A gas and A to t will be temper- 
ature dependent and hence 1/T = dScN I dEcN while 
dStot/dEtot ^ 1/T ^ dS g a S /dE gas . 

The quantity gcN should contain all the negative- 
energy bound states located in the well of the nuclear- 
plus-Coulomb potential. For positive energies, gcN 
should be independent of the bounding volume. How- 
ever apart from these constraints, there is no well jus- 
tified subdivision of g to t into its two components in the 
independent-particle model. Two methods have been uti- 
lized to calculate gcN- 
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FIG. 2: (Color online) The dependence of the compound- 
nucleus single-particle level density on the nucleon energy 
e. The displayed results have been convoluted by a Gaussian 
resolution of FHWM=150 keV. Results are shown for a) neu- 
tron in a spherically symmetric potential (Q=0) , b) neutrons 
in a deformed potential (Q=4) , and c) protons in a spherically 
symmetric potential. 



A. Subtraction Method 

In 1978 Fowler, Engelbrecht, and WoosleyQ proposed 
that g gas = gev for neutrons and thus gcN could be ob- 
tained from subtraction, i.e., gcN = gtot — gev We will 
call this the subtraction method for determining gc n and 
it has been used by m any other investigators. It is rather 
easy to show that [H El [H G3 



sub 

gcN 



(e) = $>,,(£), (17) 
id 

9td {e) = (2j + l)^( £ - £ ^) 

i 

+ -2j + l-^ 18 

7r de 

where 5ij (e) is the phase shift associated with the scat- 
tering state of energy e, orbital angular momentum £, 
and total angular momentum j. The bound-state ener- 
gies are . From Levinson's theorem[rl|. 



gij (e) de = 0. 



(19) 



It is clear that for e > 0, gij must have a net nega- 
tive contribution to balance out the positive contribu- 
tions from the bound states. However, for the large I 
waves, this negative contribution occurs at very large e 
values which are not populated in the CNE2- Thus the 
negative contributions are only important for the lowest 
I- waves. 

Near a resonance 



dSpj 
de 



T R /2 



(e - e R Y + (r*/2)' 



(20) 



where e R is the resonance energy and T R is its width. 
In the limit as — ► 0, ddej/ds — > ir6 (e — sr) and the 
resonance becomes equivalent to a bound state. There- 
fore at positive energies, g^ consists of series of res- 
onance peaks which, for low £-waves, sit on a negative 
background. 

For protons, the single-particle level density is calcu- 
lated from the nuclear phase shift, not the total phase 
shift. Hence, the subtracted level density is not actually 
g ev , the contribution from the bounding volume with- 
out any mean-field potential, as used for neutrons. In 
this case, the subtracted level density is that from the 
bounding volume containing a point-source Coulomb po- 
tential. For deformed systems, the e >0 contribution can 
be generalized as 



mcd systcn 



2ni 



-Tr 



S-\e)-S(e: 



(21) 



where S(e) is the S-matrix for scattering at energy e. 

Examples of g^V^ (e) calculated for 160 Yb are shown in 
Fig. The negative background is clearly observable in 
Fig. for neutrons in a spherically-symmetric potential. 
Here g^ (e) is negative between the resonance peaks 
for e <6 MeV. In deforming the potential, degenerate 
resonances are split and this often leads to a filling up of 
the negative background so that g^ seldom drops below 
zero. For the deformed example in Fig. |2Jd, g & ^ is only 
negative at e ~0. For protons (Fig. [21;), only very narrow 
resonances are observed well below the Coulomb barrier 
(e ~9 MeV). 
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B. Gamov Method 

By placing the CN in a bounding box, we have pro- 
duced an equilibrium model. Now within this equilibrium 
model, any arbitrary subdivision of gtot into gcN and 
g gas can be considered. Whatever subdivision is made, 
the inverse cross section Oi nv must be chosen to describe 
the absorption of nucleons from the "gas" phase space 
into the "compound-nucleus" phase space. However, in 
order for the equilibrium evaporation rate to be equated 
to that of an isolated CN which is a nonequilibrium prob- 
lem, gcN should be chosen such that f(e)gcN (e) also de- 
scribes the nucleon energy-density when there is no gas 
present. As g^ can be negative for some energies, the 
number of nucleons ascribed to the CN in these single- 
particle levels is also negative. In the equilibrium model, 
this does not pose a problem as we have an accompany- 
ing gas in the box. The total number of nucleons in any 
energy range is always positive, i.e., gtot >0. It is just our 
prescription of dividing up g to t into gcN and g gas that 
gives rise to this problem. However in the absence of the 
gas, it is not clear what physical significance a negative 
value of gcN has. 

Another criticism of the subtraction method is that it 
includes both narrow (long-lived) and wide (short-lived) 
resonances. It has been suggested that only resonances 
of lifetime longer than, or comparable to, the compound- 
nucleus lifetime should be included \PA lisl ] as the occu- 
pancy of the shorter-lived levels will not be maintained 
before the CN decays. 

Consider the analytical continuation of the S-matrix 
into the complex-energy plane. Poles of the S-matrix at 
e—e R — iT R /2 correspond to exponentially decayin g so - 
lutions to the Schrodinger equation (if T R >0)\13L Il9|. 
These are also known as Gamov |2f| or Siegert states. 
For those poles close to the real axis on the unphys- 
ical sheet, these states have a close association with 
resonances. Weidenmiiller [l7j suggested the compound- 
nucleus single-particle level density for e >0 should be 
the density of sharp resonances or Gamov states. One 
should therefore introduce a cutoff or maximum width 
Tq of the Gamov states that contribute to gcN- Now, if 
the single-particle potential is modified, for example by 
deformation, the location of the poles will move in the 
complex-energy plane. Some will become narrower and 
some wider and a number of these will cross the cutoff 
region causing discontinuities in the level density as it 
evolves with deformation. In order to avoid such dis- 
continuities, a smooth exponential cutoff of the Gamov 
states was implemented; 



ef ) exp 



-r 



(22) 



The summation includes all poles associated with bound 
states (T R =0) and Gamov states (T R >0). With this def- 
inition, g C N is always positive and thus avoids the ambi- 
guities associated with negative values as in g^}%- In the 



following sections, both definitions of gcN will be used 
to see how they affect the excitation-energy, deformation, 
and asymmetry dependencies of the level density. 



IV. COUPLED-CHANNELS CALCULATION OF 
SINGLE-PARTICLE LEVELS 

A. Theory 

In order to calculate the single-particle level densities, 
the Schrodinger equation must be solved to determine the 
bound, Gamov, and scattering states. Protons and neu- 
trons are assumed to move in an axially-symmetric mean- 
field potential which is the sum of the nuclear, Coulomb, 
and spin-obit components, i.e., 



V{r)=V N (r) + V c {r) + V so (r). 



(23) 



The nuclear potential is expressed in terms of the Fermi 
function f(x) — [1 + exp(a;)] as 



V N (r) 



■R(Q,Q) 
d(6) 



(24) 



where R (Q, 9) defines a spheroidal surface with the same 
volume as a sphere of radius R^°\ The deformation is 
express in terms of the relative quadrupole moment Q 
related to the radii, rii and r i , perpendicular and parallel 
to the symmetry axis by 



m 



Q 



Sn 



(i - 
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[R(°)y 



(25) 



The quadrupole moment is positive for prolate shapes, 
negative for oblate, and zero at sphericity. As a cali- 
bration point, Q=3.2 corresponds to a "superdeformed" 
prolate shape with the length of the major and minor 
axes differing by a factor of 2. 

The diffuseness of the nuclear potential is assumed to 
be constant perpendicular to this surface, i.e., 



d(0) 




(26) 



The deformed spin-orbit interaction can be expressed in 
terms of the momentum p and spin s operators as|22j 



V so (r) = 



V/(- 



— R so {Q, 0) . 
d so (O) ' 



x p • s (27) 



where d so is defined in terms of R so in an equivalent 
manner as in Eq. 1261 The Coulomb potential is approxi- 
mated as that from a sharp-surfaced spheroid of equiva- 
lent spherical radius Rc using the analytical expressions 
of Refs. HHH. The parameters V^ 0) , V s ( o 0) , Rf-°\ Rf \ 
R c , and are taken from the "universal" parametriza- 
tion of Ref. $2j%. 
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The solutions to the Schrodinger equation H^> = 
are obtained by expressing the wavefunction as sums of 
spherical waves specified by \i\jm). Here I is the orbital 
angular momentum, j is the total angular momentum 
and m is its projection on the symmetry axis. This latter 
quantity is conserved in the axially-symmetric potential. 
Thus 



where Uj>j m (r) are the radial wavefunctions. After pro- 
jecting on the state \£'^j'm), the Schrodinger equation 
can be written in terms of the coupled-channels equation 



(r) = 2^ r l f 2 JW ^ 



(28) 



dr 2 r 2 



(r) + ( W W>f M + D W'f M i) u *'J'm (r) 



= 0. 



(29) 



r 



Here fc = y/2[ie/h is the wave number and is the re- 
duced mass. The matrices W and D are determined 
from the matrix elements of the interaction taken be- 
tween states specified by £,j and £',j'. The matrix 
W has contributions from all three potentials (nuclear, 
Coulomb, and spin-orbit) while D is determined only 
from the spin-orbit potential. 

The boundary conditions at the origin are tt^ m (0)=0. 
If one considers N channels and chooses N initial sets 
of the derivatives duij m /dr(0) appropriately, then after 
integrating out from the origin, one can obtain TV inde- 
pendent solutions to the coupled-channels equation. Let 
these be represented by the N columns of the N x N ma- 
trix U m {r). In matrix form, the Schrodinger equation is 
then 



ar^ dr 



U m {r) =0 



(30) 



where 



Am _ 

A e.ji>3> - 



,2 + 



k z - 



6 u ,5 jr +W% e , r . (31) 



The equation is integrated out to a radius r matc h where 
V— >0 for neutrons or, for protons, only a point-source 
Coulomb term is present. At r matc hj the solutions are 
matched to specific solutions of the Schrodinger equa- 
tions pe(r). For bound states, the matching solution must 
vanish as r^oo and thus pe(r) = y/2\k\r/iv K e+ i(\k\ r) 
or pi (r)—W_ r)e+ 1 (2 |fc| r) for neutrons and protons, re- 
spectively. Here K t+ t are the modified Bessel func- 



are the Whit- 



tions of the second kind and W 
taker functions and rj = (Z — l)e 2 /(h 2 \k\). For 
Gamov states, the matching functions are outgoing 
waves; pt,{r)=kr \je(kr) + iye(kr)] or pt(r)=Fi(ri,kr) 
—iGi(r],kr) for neutrons and protons respectively. Here 
je and yi are the regular and irregular spherical Bessel 
functions and Fg and Gg are regular and irregular 
Coulomb wavefunctions. If the calculated solution is to 



represent a bound or Gamov state, then one must be 
able match the logarithmic derivatives of ui 3m and pi at 
f = f match f° r ah channels. Any linear combination of 
the column vectors of U m can be used to achieve this 
match and it is only possible when[2(| 



dU 



til (u-y'-f(pr 



where the matrix P is defined as 



= 



Pljt'j'{r) = 5 U '5 jf p e (r). 



(32) 



(33) 



l /dr{U m y L 
satisfies the 



is called the log- 
following Ricatti 



Y 2 



DY = 0. 



(34) 



The matrix Y = dU m ■ 
derivative matrix and 
equation 

— + A 

dr 

Rather than solving the matrix Schrodinger equation 
(Eq. |3UJ|, this equation can be solved directly using the 
techniques of Refs. [27Ll2^ . l29| . In fact, it is advantageous 
to solve the Ricatti equation instead of the Schrodinger 
equation as the latter suffers from numerical instabilities 
when integrating over classically forbidden regions. 

To obtain scattering solutions, the wavefunctions must 
be matched to a combination of ingoing and outgoing 
waves at r = r ma t c h- The scattering matrix can also be 
obtained directly from the log-derivative [2?]. Defining 
the matrix elements 

Jtj£'j'{ r ) = Su>Sjji krji(kr) for neutrons 

= Sig/Sjji Fg{kr) for protons (35) 

Nij£ij/(r) — SuiSjjr kryiikr) for neutrons 

= —Sa/Sjj' Gg(kr) for protons, (36) 

the K matrix is determined by 



K 



^(rmatch)-2V(r matc h) - -j^N(r match ) 

y (r match) J (r match) - ~T J ( r match) 

dr 



(37) 
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The S matrix is derived in terms of the identity matrix 
J as 

S= (I + iKy 1 (I-iK). (38) 

The calculation of g^fy from Eq.[2]can be problematic 
near very narrow resonances. However, the level density 
convoluted with a small dispersion is more easily deter- 
mined. If F(e) is the convolution function, the convo- 
luted level density is 

/>oo 

g c (e) = / g c (e')F(e - e')de'. (39) 
Jo 

Following Sandulescu et a/. |3(j, the integral along the 
real axis can be replaced by a contour integral in the 
complex-energy plane which avoids narrow resonances. 
The contour C is chosen to follow the real axis except 
near resonances with T < 50 keV where it follows a semi- 
circular path of radius 0.2 MeV around each resonance. 
From Cauchy's theorem, the final level density is 

g c (s) =J2F(e- e£) + / g c (e')F(e - e')de' (40) 

where here e n = — is the complex energy of the 

nth avoided resonance. The convolution function F was 
taken as Gaussian with FWHM=150 MeV. This small 
resolution does not have any significant affect on the de- 
duced level densities in this work. 




FIG. 3: (Color online) Evolution of the real energy 
e R for bound (e R <0, T R =0) and narrow Gamov 
(e R > 0, F R < 1 MeV) states as a function of deformation. 
The width of the curves reflects the value of F R . Results are 
only shown for a) neutron and b) proton 771^=1/2^ states in 
190 Yb. 



B. Results 



An example of the evolution of bound single-particle 
levels and the real part of narrow Gamov states 
(r <0.5 MeV) with deformation is shown in Fig. [3] for 



m 



■i neutrons in 190 Yb. The results were obtained 



by including all channels with I <20 in the coupled- 
channels calculations. The bound states and resonances 
levels both move around with deformation, but levels of 
the same m n values avoid crossing each other. Bound 
levels that pass through e=0 immediately become nar- 
row resonances and vice versa. As the resonance energy 
increases, the width of a resonance generally increases as 
shown in Fig. 

The behavior that bound states turn into narrow reso- 
nances is quite general except if there is no barrier (cen- 
trifugal or Coulomb) .13]. For j=\ (^=0) neutrons there 
is no barrier and bound states passing through e=0 turn 
into virtual states 01 . A virtual state is associated with a 
pole of the S matrix on the real e axis at energy e v = —e v 
(e v > and small). Both bound and virtual states have 
real negative energies and purely imaginary wavenum- 
bers k. However for bound states, the imaginary part of 
k is positive while it is negative for virtual states. In fact 
when there is no barrier, both bound and virtual states 
with small energies have important influences on the scat- 
tering at small positive values of e. This has implications 



for the single-particle level density calculated with the 
subtraction method. Consider a spherically symmetric 
potential. The £—0 contribution to the single-particle 
level density from a virtual state for small e values is[l3| 



sub 

9o 




e+e„ 



if e > 0, 
if e < 0. 



On the other hand for a bound state at 
and small), the contribution is[l3| 



(41) 



-e b (e b > 



sub 

9o 



i 



5{e- 



s+e b 

eb) 



if e > 
if e < 0. 



(42) 



As e — > +0, <7Q Ub — > +00 for the virtual state and -co 
for the bound state. In the limit that e v — > and 
e b — > 0, then in both cases go ub (s) — > S(e)/2 and this 
represents half a level. Thus in the deformation region 
over which a bound level becomes a virtual state, g^V^ 
evolves smoothly. 

Only the m=i + states contain any \t—Qi) compo- 
nent in their wavefunction. For these states, the behav- 
ior as a bound state passes through e—0 is more com- 
plex. Sometimes they become narrow resonances, some- 
times they become virtual states, and other times they 
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FIG. 4: (Color online) Variation with nucleon energy e of 
the single-particle level density determined by the subtraction 
method for m lr = l/2 + neutrons in 190 Yb. Results are shown 
for two neighboring deformations Q=1.75 and 1.625. 



progress in complicated manner to a wider resonance. As 
an example, Figure 0] shows the m=i + contribution to 
the single particle level-density glj 2 + ( e ) at sman positive 
energies obtained for two neighboring values of Q. For 
Q=1.625, g{f 2+ -> -oo as e -> +0 as in Eq.gl On the 
other hand, at Q=1.75, ~~ * as m Ea. 1411 In 
this case the behavior of <7f /2+ ^ s consistent with a bound 
state at Q= 1.625 passing through £=0 and becoming a 
virtual state at Q=1.75. 

For the Gamov method, <?Jy 2 + ( e ) does n °t evolve 
smoothly when a bound state becomes a virtual state; 
bound states are always counted as a full level while vir- 
tual states are not counted at all in Eq. [221 However, 
virtual states can cause long time delays in scattering 
like narrow resonances and some thought should be 
given to expanding the definition of g^ N to include some 
contribution from virtual states and so make the evolu- 
tion with Q smoother. In the present work this is not 
a significant issue as, in the range of deformation inves- 
tigated (—2.5 < Q < 4), there is typically only 2 small 
discontinuities. 

To visualized the gross differences between the subtrac- 
tion and Gamov methods it is useful to smooth the single- 
particle level density gcN- Figure displays smoothed 
neutron and proton single-particle level densities gcN for 
160 Yb. The Strutinsky smoothing discussed in Sec. I VII 
was utilized. For neutrons, gcN peaks near e=0. The 
peak is lower in magnitude for the subtraction method 
due to the presence of the negative background. As res- 
onances at larger e R tend to have larger widths, the 
Gamov method, which excludes these wide resonances, 
makes gcM drop quickly to zero for e ^>0. The effect is 
more pronounced the smaller the cutoff width Tq. Pro- 
tons exhibit similar behavior except they peak closer to 
the Coulomb barrier whose magnitude is indicated by the 
arrow in Fig. [SJd. 
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FIG. 5: (Color online) Smoothed single-particle level den- 
sities calculated for a) neutrons and b) protons in 160 Yb. 
The curves labeled "sub" were obtained from the subtraction 
method. The other curves were obtained from the Gamov 
method with the indicated values of To in MeV. The Coulomb 
barrier for protons is indicated by the arrow in b). 



There are two important results to highlight. First 
for protons, g(e) is almost independent of the method 
of calculation for small positive energies well below the 
Coulomb barrier. These positive-energy states have the 
most influence of the level density and thus the contin- 
uum corrections for protons will generally be less impor- 
tant that those for neutrons. Secondly for neutrons, the 
inclusion of the negative background and the wide reso- 
nances in the subtraction method have opposite effects 
and partially cancel. In the end, the level densities cal- 
culated with both methods are found to be similar (see 
Sec. IV 11 At . 

The subtraction method does not give a strong defor- 
mation dependence of <?cw near sphericity. As an exam- 
ple, the smoothed single-particle level density for neu- 
trons in 160 Yb is plotted in Fig. & for Q=-0.5, 0, 0.5. 
The curves for all three deformations lie almost on top 
of each other. In contrast, the Gamov method exhibits 
a strong dependence. The results, displayed in Fig. 03 
for ro=l MeV, show that gcN decreases in magnitude 
for e > as a deformation is imposed on the compound 
nucleus (either prolate or oblate). To investigate this 
behavior, let us concentrate on the splitting of Gamov 
states which are degenerate at sphericity. For example 
in Fig. the evolution with deformation of the energy 
and width for a group of Gamov states associated with 



9 



- 6 
> 



Z 
'60 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1"| 

a) Subtraction Method 
Q 

0.0 

0.5 

-0.5 




-50 -40 -30 -20 -10 10 20 30 
e (MeV) 

FIG. 6: (Color online) Smoothed single-particle level densi- 
ties calculated for neutrons in 160 Yb at three defomations; 
Q=— 0.5,0, and 0.5. In a) results obtained with the subtrac- 
tion method are displayed while in b), the Gamov method 
with Fo=l MeV was used. 




FIG. 7: (Color online) Curves show the variation with defor- 
mation Q of the energy e R and width Y R of the Gamov states 
associated with the j=19/2 _ neutrons in 160 Yb. The solid 
data points indicate the average energy and width for these 
states. 



j=19/2~ neutrons in 160 Yb is displayed. In Fig. 03, the 
real energies e R of the Gamov states fanout with increas- 
ing deformation (both oblate and prolate). On average, 
the mean value of e R changes very little with deforma- 
tion. These mean values are indicated by the data points. 
In contrast, the widths T R of Gamov states show a differ- 
ent behavior in Fig. 03. Although a few of the states for 
oblate deformations show a reduced width compared to 
sphericity, by in large the widths of most states increase 
with deformation. The average widths, indicated by the 
data points, have a minimum at sphericity. This behavior 
is typical of the splitting of all degenerate Gamov states 
and therefore the average increase in these widths with 
deformation reduces their contribution to g^ N (Eci.l22jl. 
This explains the observed deformation dependence dis- 
played in Fig. The strong dependence of g^ N on 
deformation is reflected in the deformation dependence 
of the level density (see Sec. IVI|) . 



V. NUCLEAR LEVEL DENSITY WITH 
PAIRING 

The simple discussion of the level density in Sec.lTTlfor 
a single particle type and no interactions is extended in 
this section to include both neutrons and protons and 
the pairing interaction. The grand potential of a two- 
component Fermi gas is the sum of the proton and neu- 
tron contributions, i.e., 



SI (a n , a p , 13) = Q n (a n ,/3) + fl p (a p , 13) 



(43) 



The pair ing interaction is considered in the BCS 
modelpll l32l l33l |. In this model, the grand potential 
is related to the grand partition function (f2„ (a„ , (3) — 
— \nZ n / (3) and for neutrons it is given by 

fi» («»,/?) = ^+ J 9 -4^-[e-» n -E]de 

_ 2 fg n (e) 

0, 



■ In [1 + exp {-(3E)\ de (44) 



where (3 = 1/T, T is the temperature, and /i n = a n / (3 is 
the chemical potential. The quasiparticle energies are 



E 



(45) 



The gap parameter A n is determined from the gap equa- 
tion 



2 

G n 



g n (e) tanh (J,) 
2 E e 



(46) 



where G n is the pairing strength. The level density at 
an energy E ' can be obtained from the inverse Laplace 
transform of the grand partition function 
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^ p-\-too n-f-lOO r-f-tOO 

E tot ,N,Z) = — / / / Zexp(-a n N - a p Z + (3E) da n da p d/3 

(27Tz) J — too J — zoo J — too 



(47) 



which can be evaluated approximately by the saddle- 
point method to give 

p(Et°\N,Z)= e ^[ S) (48) 

Here the energy (E tot = + £* ot ), entropy (S = S n + 
S p ), and particle number are determined from the follow- 
ing equations: 



_g»(e) 



e — 



Gr, 



S n — 



3n(e)ln 
9n{e 



E 



( E 
1 + exp - — 



ds 



(49) 



T 



E 
T 



N 



9n(e) 



l + exp(f) 



■<fe, 



g - 
E 



tanh [ — 



ds 



2T J 



(50) 
de. (51) 



The quantities E p ot , S p , and A p for protons are obtained 
from similar expressions and the determinant D is now 



D 



tr InZ d 1 \nZ 
d 2 InZ 



da n da p 
d 2 InZ 



da n da p 
d 2 \nZ 

da* 
d 2 \nZ 



d 2 In Z 
da n d/3 
d 2 InZ 
da p d(3 
d 2 InZ 



da n d/3 da p df3 dp 1 



(52) 



Expressions for the evaluation of this determinant in 
terms of the single-particle level densities can be found in 
Ref. [32j • At some critical temperature T crlt , the gap pa- 
rameter vanishes and the excitation energy and entropy 
are those of a noninteracting Fermi gas, i.e., Eqs. [21 0] 
andH 



VI. DEFORMATION ENERGY 

The level density will be calculated as a function of 
excitation energy. The excitation energy is given in terms 
of the thermal contribution E th =E tot (T, Q) - E tot (0, Q) 
and the deformation energy Ed e f{Q)- In the Strutinsky 
procedure [34J, the deformation energy is given by two 
terms 



EdeAQ) - SE(Q) + V def (Q) 



(53) 



where the liquid-drop deformation energy Vdef{Q) de- 
scribes the average deformation energy with shell oscilla- 
tions averaged out. The corrections 5E(Q) to the liquid- 
drop energy are determined from the single-particle levels 



and have contributions from both neutrons and protons, 
i.e., SE = SE n + 5E P . Following Ref. H^, we define the 
shell corrections as 



6E k (Q)=E% 



f (0,Q) 



Et ot (o,Q) 



(54) 



where k=n or p, and Ef?* is the total energy determined 



with pairing (Eq. |4*9"|) . 
particle level densities 



but with the smoothed single- 



g(e) = / g(e')F(e - e')de' . 



The smoothing function used is 



F{e) 



exp 



c P ( e -) 



(55) 



(56) 



where the smoothing range 7 must be taken to be of the 
order of the intershell separation in order to washout the 
oscillations. The curvature correction of order p~2M is 



A I 



C p (x) = J2 



n=0 



(-1)" 
2 2n n \ 



H 2n (x) = 4f (x 2 ) . (57) 



This curvature correction is included to provide selfcon- 
sistancy for g(e), i.e., a smoothed function should not be 
affected by the smoothing procedure. Thus if g(e) is a 
polynomial of order 2M+1 or lower, it will be unchanged 

1 /2 

after the smoothing. The functions H n and L 2p are Her- 
mite polynomials and associated Laguerre polynomials, 
respectively. 

In the original Strutinsky smoothing procedure, the 
smoothing parameter s 7 and p are chosen to satisfy the 
plateau condition |34t l35| 



d~, 



dp 



0, 



(58) 



over some range in both 7 and p. Thus in this range, the 
shell correction should depend neither on the smooth- 
ing range or the order of the curvature correction. The 
plateau condition can be satisfied for single-particle lev- 
els associated with infinite potentials such as a harmonic 
oscillator or an infinite square well. However for a finite- 
depth potential, such as those considered in this work, 
the plateau condition is often not met, i.e., one cannot 
find a re gion of 7 and p over which the shell correction is 
constant |36l 1371 138| . 

An alternative procedure from Refs. H^H^ was tried, 
however this was found to problematic in some cases. In- 
stead the method that is used in this work relies on the 
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FIG. 8: (Color online) Variation of the neutron shell correc- 
tion to the deformation energy determined for 170 Yb. a) The 
absolute correction obtained with the indicated smoothing 
ranges 7. b) The data from a), but now the mean correction 
over all calculated deformations {8 (Q)) have been subtracted 
out. 



observation that the relative correction for different de- 
formations is independent of 7 and p once 7 has a value 
above ~1.2 fko. The actual 7 value at which the rel- 
ative correction plateaus depends on the order p used. 
However once the plateau is reached, the relative correc- 
tions are independent of p. As an example, the correc- 
tion factors for neutrons obtained with p—12 for various 
7 values are plotted in Fig. [SJl. The absolute values of 
these corrections vary continuously with 7 and do not 
plateau. However, it does have a minimum in the interval 
2.0< 7 <2.5 in this example. Now apart from the dashed 
curve obtained with 'y—huj, all other 5 n {Q) curves have 
almost the same shape indicating the relative shell cor- 
rection is constant. To highlight this, the average value of 
the correction over all deformations (6 n (Q)) is subtracted 
out for each smoothing range. The remaining correction 
S n {Q) — (6 n (Q)) is plotted in the Fig.[Sj3. All the curves 
for 7 > fko now collapse to essentially a single curve. 

If we cover a large enough range of deformations, the 
average shell correction (Sk(Q)) is expected to be zero, 
thus we have taken the values plotted in Fig. [SJ} to be 
absolute corrections. Thus the appropriate smoothing 
range is the value which causes (5k(Q))—0. We expect 
the range of deformations explored in this work (-2.5< 
Q <4.0) to be adequate as we always see at least a couple 
of "oscillations" in 5k and thus expect the average to 
reflect the true average over all deformations. 



FIG. 9: (Color online) Variation of four calculated energies 
with deformation for 170 Yb and 150 Yb. The curve E def is the 
deformation energy with shell corrections. For comparison 
the liquid-drop deformation energy Vu is also shown. The 
higher two curves give the energies corresponding to the crit- 
ical temperatures T crtt for protons and neutrons. The thick- 
dashed lines indicate the deformation which maximized the 
level density (subtraction method) for each excitation energy. 



Finally the excitation energy is measured with respect 
to the ground-state energy, i.e., the total excitation en- 
ergy is 

E* =E th + E def (Q)-5W. (59) 

Here the shell correction 5W represents the difference 
between the liquid-drop and the minimum or ground- 
state deformation energy, i.e., SW = min [Ed e f(Q)]- Note 
that any error in the absolute value of 5k(Q) affects both 
Edef(Q) and SW equally, and therefore the excitation 
energy is not sensitive to the absolute shell correction. 

The liquid-drop deformation energy is taken 
from Refs. |3^, |4jj. The gap strength G„ iP for 
neutrons and protons is determined from setting 
A^p(Q=0,T=0)=12/VA MeV 6]. Here is the gap 
parameter obtained from Eq. 1461 with the smoothed 
single-particle level densities g^p. Examples of the 
deformation energy are shown in Fig. [5] for systems with 
deformed ( 170 Yb) and spherical ( 150 Yb) ground states. 
Also shown are the excitation energies corresponding to 
the critical temperature T^™* for neutrons and protons 
where the pairing gap vanishes. 
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TABLE I: Nuclei studied in this work and the value of the 
shell correction 5W , minimum cost E™'!,\, and pairing factor 
5P. 



Nucleus 


SW 


jpmm 

^cost 


SP 




(MeV) 


(MeV) 


(MeV) 


4,1 Si 


1.1 


1.9 


2.7 


40 g 


0.6 


5.4 


2.7 


40 Ar 


0.1 


9.4 


2.7 


40 Ca 


-0.7 


7.4 


2.7 


40 Ti 


0.6 


3.1 


2.6 


60 Ti 


-0.8 


3.1 


2.6 


60 Cr 


0.2 


5.7 


2.6 


60 Fe 


0.4 


8.4 


2.5 


60 N j 


-0.6 


10.1 


2.5 


60 Zn 


-0.6 


6.9 


2.5 


60 Ge 


-1.6 


3.8 


2.6 


170 Nd 


0.6 


2.5 


2.4 


190 Yb 


-1.6 


3.8 


2.4 


180 Yb 


0.5 


5.6 


2.3 


170 Yb 


-0.8 


7.6 


2.3 


160 Yb 


0.7 


9.9 


2.2 


150 Yb 


-2.4 


8.9 


2.2 


170 Hg 


-2.9 


10.0 


2.2 


238 tj 


-0.7 


5.5 


2.3 



VII. CALCULATIONS 

A. Excitation-Energy and n-p Asymmetry 
Dependencies 

The level density as function of excitation energy was 
calculated for even-even nuclei with A=4Q and 60. All 
such nuclei with E™£ t >1.9 MeV were included in the 
study. These include 40 Ti and 60 Ge which are just be- 
yond the proton-drip line. Calculations were also per- 
formed for five even-even Yb nuclei from 150 Yb to 190 Yb 
covering the range of n-p asymmetry centered on /3-stablc 
nucleus 170 Yb. In addition two other A—\7Q nuclei, 
170 Nd and 170 Hg, with extreme values of n-p asymmetry 
were included. Again, 150 Yb and 170 Hg are just beyond 
the proton drip line. Finally the heavier /3-stable system 
238 U was also studied. All nuclide studied are listed in 
Table U along with their values of SW and -E™ s ". 

At each deformation Q, the level density and excita- 
tion energy are calculated for an array of temperatures 
each separated by 0.05 MeV. The level density for a given 
excitation energy is then obtained from interpolating be- 
tween there results. Subsequently, the deformation of 
the nucleus at each excitation energy is determined as 
the value which maximizes the level density. As an ex- 
ample, the deformation as function of excitation energy is 
plotted in Fig. as the thick-dashed curves. For the de- 
formed ground-state system 170 Yb, the deformation de- 
creases with excitation energy and vanishes at E*=70 
MeV. The spherical ground-state system 150 Yb remains 
spherical at all excitation energies. 

The variation of the resulting level density with exci- 
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FIG. 10: (Color online) Variation of the nuclear level density 
with excitation energy for a) the seven A ~170 nuclei, b) 
the six A— 60 nuclei, and c) the five A=40 nuclei. All are 
calculated up to the excitation energy where T=E™f t . These 
values are indicated by the solid circular symbols if they are 
within the displayed ranges and are labelled. 



tation energy obtained with the subtraction method is 
plotted in Figs. HUfa . b, and c for A ~170, A=60 and 
A=40, respectively. In Fig. HOfa , where A is not con- 
stant, the quantity log (p^4 5 / 4 ) A 1 / 2 rather than log (p) 
has been plotted to account for the A dependence based 
on the Fermi-gas formula with a oc A. Curves for all nu- 
clei are only extended up to the excitation energy where 
T=£'™ I S ". It is clear from this figure, that the level den- 
sity has no substantial dependence on n-p asymmetry, all 
curves with similar A values practically overlap. Similar 
conclusions were also obtained with the Gamov method. 
For example, the level densities for A ^170 and A=4Q 
are shown in Figs. ITT| and 1121 respectively, for ro=1.0 
and ro=0.01 MeV. Again, the curves for similar A val- 
ues fall almost on top of each other. 

For T > T^"'and T > T^ rit in even-even nuclei, the 
excitation energy is often backshifted by the conden- 
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FIG. 11: (Color online) Same as for Fig. HOh ,. but now the 
level density for A ~ 170 is calculated with the Gamov method 
using the cutoff decay widths of ro=l and 0.01 MeV. Only 
curves that are distinguishable from the others are labelled. 



FIG. 12: (Color online) Same as for Fig. HUfc . but now the 
level density for A=40 is calculated with the Gamov method 
using the cutoff decay widths of ro=l and 0.01 MeV. Only 
curves that are distinguishable from the others are labelled. 



sation energy when comparing level densities |4l| . i.e., 
U=E* — SP where the smoothed condensation energy is 

1 ~2_ 1 ~2_ 

SP = g A « AViOn) + 2 A p dpiVp)- (60) 

For odd-even and odd-odd nuclei, SP should include the 
pairing correction in the semicmpirical mass formula. In 
the comparison of level densities in Figs. II OH 121 the role 
of pairing is not important as the condensation energy is 
relatively constant for each mass region (see Table HJl. At 
high excitation energies where shell effects are expected 
to be washed out, the excitation energy is also shifted 
by the shell correction Thus at high energies, 

a shifted Fermi-gas expression is often assumed. In this 
case, the entropy is 

S = 2yJa(E* - 5P + SW) (61) 

where a is the asymptotic level-density parameter. 

To see whether this formalism is consistent with the 
calculations of this work, the asymptotic level-density 
parameter is calculated from Eq. 1611 for all excitation 
energies. Examples of the resulting level-density param- 
eters are displayed in Fig. ^5] for the subtraction method 
and in Figs. 1141 and 1151 for the Gamov method. Above 
(E* -5P + SW) I A >0.3 MeV where shell and pairing 
effects are expected to the quenched, the deduced level- 
density parameter is rather constant. Note, Figs. I13HT5I 
have offset origins on the y axis to accentuate the differ- 
ence between the different nuclei. Quite surprisingly, the 



inclusion of realistic single-particle level densities includ- 
ing continuum corrections does not cause strong devia- 
tions from the basic Fermi-gas expression which was de- 
rived for constant g (e). A similar conclusion was found in 
Hartee-Fock calculations of 208 Pb performed by Bonchc, 
Levit, and Vautherin|43j. 

The dependence of a on excitation energy is not com- 
pletely flat, all calculations show some small negative 
slope as is expected in a lowest-order expansion in tem- 
perature (Eq. 111(1 . This is most significant in the Gamov 
method (Figs. 114b and 115b ) with the smallest value of 
cutoff the width (r =0.01 MeV). The plots of the asymp- 
totic level-density parameter also highlight the small dif- 
ferences between the nuclei with the similar A values 
which were difficult to see in Figs. llOH12l The differences 
between the curves would become more significant if they 
are extended beyond the point T > E^™\ where the the 
statistical model is problematic. For A—A0 : the asymp- 
totic level-density parameter can be somewhat smaller for 
the systems with extreme n-p asymmetries. For example, 
look at the results for 40 Si in Fig. 113b with the subtrac- 
tion method and for 40 Ti in Fig.[]2b with r =0.01 MeV. 
For both of these cases, the asymptotic region is only 
approach when E™"\ » T. To investigate the small n-p 
asymmetry dependence more systematically, the value of 
a determined at (E* - SP + SW) /A=0.5 MeV is plotted 
in Fig. ^3 as a function of N — Np (A), the distance from 
the /3-valley of stability. Results are shown for the sub- 
traction method (filled circles), and the Gamov method 
with three values of the maximum width; Fq=1 MeV 
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FIG. 13: (Color online) Level-density parameters a deduced 
from the calculated entropy and Eq. 1611 as a function of the 
shifted excitation energy per nucleon. Results were obtained 
using the subtraction method to treat the continuum. The 
data points indicated on each curve are the point at which 



(hollow squares), ro=0.1 MeV (filled diamonds), and 
ro=0.01 MeV (hollow triangles). Generally, the deduced 
values of a are approximately independent of which treat- 
ment of the continuum was used and are almost constant 
for each mass region. In the Gamov method, the nuclei 
with the extreme values of N — Np (A) show the greatest 
sensitivity to Tq. In this case the values of a obtained 
with ro=0.01 MeV are often slightly smaller. 

The mass dependence of the level-density parameter 
at (E* -SP + SW) /A=0.5 MeV is displayed in Fig. [T7| 
as the data points. The extracted points were fit by the 
commonly used formula 

B s a s A^ 



a 



a, 



A 



(62) 



which includes volume and surface contributions where 
a v and a s are the coefficients for these two quan- 
tities. The dimensionless parameter B s gives the 
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FIG. 14: (Color online) Same as for Fig. Il.'-ih . but now the 
results for A ~170 were obtained from the Gamov method 
with a) To=l MeV and b) Fo=0.01 MeV. Only curves that 
are distinguishable from the others are labelled. 



surface area of the nucleus relative to its spherical 
value. As all the systems studied are spherical at 
(E* -SP + SW) /A=Q.5 MeV, B s was set to unity. The 
dashed curve in Fig. El shows the fit obtained with 
Eq. E2 The fitted coefficients are a„=0.078 MeV -1 
and cv s =0.146 MeV -1 . For comparison, curves for 
the level-density parameters from Toke and Swiatecki 
(a„=0.068 MeV -1 , a s =0.274 MeV _1 )Q and Ignatyuk 
et al. (a„=0.073 MeV" 1 , a s =0.095 MeV _1 )0| are also 
displayed. The fitted surface coefficient is intermediate 
in value between these two other prescriptions, but closer 
to Ignatyuk et al. 

Many experimental studies have adopted the 
excitation-energy dependence of the level density 
suggested by Ignatyuk et o/.0, which includes the 
washing out of shell effects with increasing temperature. 
The entropy is expressed in terms of an excitation-energy 
dependent level-density parameter, i.e., S = 2y/a(U)U 
where 



i(U) = a 



l + h(U) 



SW 
IT 



(63) 



The function h, determining the behavior at low excita- 
tion energies, is given by h(U)=l-exp (— *yU). The pa- 
rameter 7 gives the energy scale over which shell effects 
are washed out. At high excitation energies where h — ► 1, 
the Ignatyuk formalism leads to the expected dependence 
of Eq.|^ To determine how well Eq.[U3|can describe the 
entropy calculated in this work, asymptotic level-density 
parameters were determined at each excitation energy 
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FIG. 15: (Color online) Same as for Fig. 113b . but now the 
results for A— 40 were obtained from the Gamov method with 
a) r =l MeV and b) r =0.01 MeV. 
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FIG. 17: (Color online) Asymptotic level-density parameter 
a as a function of nucleon number A. The data points are the 
values determined from the calculations of this work. The 
dashed curve shows a fit to these data. Curves are also shown 
for the prescription of Toke and Swiatecki|44l| and Ignatyuk 
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from Ea. 1631 and the parameter 7 was adjusted to mini- 
mize the spread of these deduced a values at low excita- 
tion energies. For the A ~170 nuclei, 7 was determined 
by this procedure as 0.035 MeV -1 and the a values are 
plotted in Fig. 118b . These are to be compared to the 
a values obtained for 7=00 in Fig. 118b ,. The condition 
7=00 corresponds to h=0 and the latter values are the 
same as deduced from Ea. 16 Hand plotted in Fig. 113b. The 
spread in the a values at low excitation energies observed 
in Fig. I18fa are almost removed in Fig. 118b and thus this 
indicates that the Ignatyuk formalism describes the fad- 
ing out of shell effects adequately for this mass region. 
The deduced value of 7 is of similar magnitude to the 
value 0.05 obtained by Ignatyuk et al. 45] by fitting neu- 
tron resonance data. Schmidt et al.\ 16 have extracted 
a mass-dependent value of 7 and, for A— 170, they find 
7=0.045. Again of similar magnitude to the value of this 
work. 

For the lighter mass regions (A=A0 and 60), a similar 
reduction in the spread of the deduced a values was not 
achieved. Thus for these light systems, the description of 
the level density in the region where shell effects are still 
important is more complex than this Ignatyuk treatment. 



B. Deformation Dependence 



FIG. 16: (Color online) Calculated level-density parame- 
ters at (E* -SP + 5W) /A=0.5 MeV are plotted verses N - 
N/3 (A), the neutron number separation from the j3 valley of 
stability . The data points were obtained with the subtrac- 
tion and Gamov methods. For the latter, the three indicated 
value of the cutoff width Fo were used. For ^4=40 and 60, the 
fitted variation of the level-density parameter (cases B and C) 
from Ref. Q are shown by the dot-dashed and dashed curves, 
respectively. 



The level density of strongly deformed nuclei is of in- 
terest. In the statistical model, the fission decay rate is 
determined from the level density of the deformed saddle- 
point configuration. The deformation dependence of the 
level density is also needed to determine the equilibrium 
shape distribution of compound nuclei [3^|. This distribu- 
tion can be important in determining the emission rates 
of a and heavier fragments |47j. The deformation depen- 
dence for 170 Yb at various excitation energies is displayed 
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FIG. 18: (Color online) Asymptotic level-density parameters 
determined for the A ~170 systems with Eq. 1631 a) The 
parameter 7 was set to infinity. These results are the same as 
those shown in Fig |13h . b) The value 7=0.035 was obtained 
from minimizing the spread in the curves at low excitation 
energy. 



in Fig. To highlight the deformation dependence, 
the level densities are normalized to the maximum value 
for that excitation energy. The results obtained with 
the subtraction method are shown as the solid curves, 
while the dashed curves are from the calculation with 
the Gamov method (To=l MeV). At the lowest excita- 
tion energies, the level density is largest for deformations 
close to the ground-state value (Q=0.875 for 170 Yb). For 
this nucleus at these excitation energies, the continuum 
is not sampled significantly and the results for the two 
methods are almost identical. At an intermediate en- 
ergy (~100 MeV), shell effects have melted and the level 
density peaks for spherical shapes but the distribution is 
quite broad. Again, the results are similar for the two 
methods. However at higher excitation energies where 
the continuum becomes more important, the results ob- 
tained with the two methods are quite different. For 
the subtraction method (solid curves), the dependence 
on deformation near sphericity decreases. The curves 
become broader with increasing excitation energy as ex- 
pected when the temperature increases. Contrary to this, 
the dashed curves obtained with the Gamov method be- 
comes narrower. The underlying reason for this behavior 
can be traced to the variation of the resonance widths 
with deformation displayed in Fig. {7\ and discussed in 
Sec. IIV Bl As a spherical system is deformed, then on 
average, the widths of the resonances increase leading 
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FIG. 19: (Color online) Deformation dependence of the level 
density calculated for 170 Yb at the indicated excitation en- 
ergies. The level densities are normalized to the maximum 
Pmax for that excitation energy. For clarity the results for 
each successive excitation energy are shifted up along the y 
axis. The thin lines in each case correspond to the shifted 
x axis. Results obtained with the subtraction method are 
shown as the thick-solid curves, while the dashed curve were 
obtained from the Gamov method with Fq=1 MeV. 



to a decrease in the single- particle level density in the 
Gamov method. Thus, this behavior leads to a favoring 
of spherical shapes. 

The strong peaking of the level density at sphericity 
for the Gamov method is quite general. Results are 
shown in Fig. HJ for 160 Yb and 60 Ni at E*/A =1.875 
MeV. In both cases, the continuum contributions are sig- 
nificant. Curves are shown for the subtraction method 
(solid) and for r =l MeV(dashed) and r =0.01 MeV 
(dot-dashed). For both nuclei, the level density distri- 
butions obtained with the Gamov method are narrower 
than those from the subtraction method. However, the 
peaking at sphericity is even stronger for smaller values 
of the cutoff width To. Also, the effect is stronger for the 
heavier system. These dependencies are quite general. 

The degree to which the spherical shape is favored 
in the Gamov method also depends on the n-p asym- 
metry. Figure 1211 displays the deformation depen- 
dence of the level density for three Yb isotopes all at 
E* /A—1A8 MeV. The dashed curves, calculated with 
the Gamov method, indicate that the effect is stronger 
for the very neutron-rich 190 Yb isotope. Because of the 
small neutron separation energy for this system, the im- 
portance of the positive-energy neutron levels is much 
greater than for the /3-stable 170 Yb system. The proton- 
rich 150 Yb system shows an even smaller effects than the 
/3-stable nucleus. The most important positive-energy 
proton levels are below Coulomb the barrier and thus 
they are all narrow compared to the cutoff width Tq . As 
a general rule, modifications induced by deformation are 



17 




FIG. 20: (Color online) Deformation dependence of the nor- 
malized level density as in Fig. EH Results are shown for 
160 Yb (£*=300 MeV) and 60 Ni (£*=112 MeV) with the sub- 
traction method (solid curves) and with the Gamov method 
for r = l MeV (dashed curves) and 0.01 MeV (dot-dashed 
curves) . 



of lesser importance for protons compared to neutrons. 
Of course for very light nuclei, the Coulomb barrier is 
smaller and its ability to suppress these deformation ef- 
fects is reduced. 

The deformation dependence from the Gamov method 
has important consequences for fission. The fission decay 
width is determined from the ratio of the level densities 
at the saddle-point and equilibrium configurations. The 
favoring of spherical nuclei at high excitation energies 
will lead to a suppression of the fission width and thus 
an increase in the probability of the competing evapo- 
rative decay modes. The total fission cross section will 
therefore be reduced and, for events that do fission, it 
will occur later in the decay cascade. Such effects have 
been observed experimentally. Measurements of pre and 
postscission multiplicities of evaporated particles in co- 
incidence with fission fragments have shown that these 
fragments are created with little excitation energy even 
when the initial CN excitation energy is large 48] . The 
standard interpretation of these results is in term of dy- 
namical effects (43, but it is clear that if one adopts the 
Gamov method, then it can explain part, or possibly 
most, of the experimental observations. The magnitude 
of the predicted effect will depend on the value of To. 

Fission is most important for the heavier systems and 
thus it is of interest to examine the deformation depen- 
dence determined for 238 U. This in shown in Fig. [23 
for three excitation energies. For the lowest value 
(E* = 16 MeV), the level density is again largest for defor- 
mations around the ground-state value. There is again 
no difference between the subtraction (solid curve) and 
Gamov (dashed curve) methods. (The curves are indis- 




FIG. 21: (Color online) Deformation dependence of the nor- 
malized level density for three Yb isotopes at the indicated 
excitation energies. Solid and dashed curves were obtained 
with the subtraction and Gamov (ro=l MeV) methods, re- 
spectively. 



tinguishable) . At i?* = 140 MeV shell effects have melted 
and, in both methods, the level density is largest for 
spherical systems. Only a small difference between the 
two methods is observed. Again at the highest excitation 
energy (£T=560 MeV), the two methods give very differ- 
ence results. The Gamov method is strongly peaked at 
sphericity. In contrast now, the results with the subtrac- 
tion method show this nucleus is unstable with respect to 
prolate deformations, i.e., the level density increases with 
increasing Q. This is basically a fission instability, how- 
ever to fully treat fission one should include more shape 
degrees of freedom. If the level-density parameter is de- 
formation dependent as in Eg. 1621 then the fission barrier 
is temperature dependent |49|. The result obtained with 
the subtraction method therefore represents the situation 
where the temperature-dependent fission barrier has van- 
ished. In terms of level density, there is no saddle-point, 
i.e., a configuration of low level density which represents 
a bottleneck through which the system must pass in or- 
der to fission. Therefore, fission stability can be quite 
different for the two methods of treating the continuum. 

As the deformation-dependence of the level-density pa- 
rameter plays an important role in fission, the applicabil- 
ity of Ea.l62lwas investigated for the subtraction method. 
(It is clearly not applicable for the Gamov method). For 
a given excitation energy and deformation, the asymp- 
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FIG. 22: (Color online) Deformation dependence of the nor- 
malized level density for 238 U as in Fig. 1191 
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FIG. 23: (Color online) Asymmptotic level-density paramters 
deduced for deformed 150 Yb nuclei (oblate and prolate 
shapes) as a function of the relative surface area B 3 . Results 
are shown for the indicated thermal excitation energies. For 
clarity, results at some excitation energies have been shifted 
up along the y axis by the indicated amounts. 



totic level-density parameter was determined as 



S = 2Ja (E* - SP + SW - Vdef(Q)) 



(64) 



assuming shell and pairing effects are washed out. This 
differs from Eq. 1611 in that now the liquid-drop defor- 
mation energy has also been subtracted from the ex- 
citation energy. Results obtained for 150 Yb are dis- 
played in Fig. 123 as the data points for various values 
of E t h — E* — SP + SW — V de f(Q), the asymptotic ther- 
mal excitation energy. The extracted values are plotted 
against B s , the relative surface area (B s =l is sphericity). 
Except for the lowest value of E t h where shell oscillations 
are still present, they increase almost linearly with B s . 
The solid curves display linear fits to the extracted val- 
ues and, from the fitted slopes, the surface coefficient a s 
(Eg. I62[) can be determined. Ignoring E t h=70 MeV, the 
slopes, and thus the a s coefficients, are almost indepen- 
dent of excitation energy. The a s coefficients obtained 
from all calculated A ^170 nuclei are plotted in Fig. 1241 
Apart from the lowest excitation energies where shell ef- 
fects are still important, all the a s coefficients are similar, 
almost independent of excitation energy and n-p asym- 
metry. The values of a s determined by this procedure 
are quite similar to the value obtained from fitting the A 
dependence of a in Sec. lVII A"l . This value is indicated by 
the dotted line in Fig. [21 Clearly Eg. 1621 provides a rea- 
sonably consistent description of the A and deformation 
dependencies of the level-density parameter for the sub- 
traction method. The small difference between the ex- 
tracted values of a s from the two procedures may be due 
to the fact that curvature and higher-order corrections to 
the level-density parameter have been ignored. Thus, 
these must be small (at least for the deformations con- 
sidered) in order to get such good agreement from the 
two procedures. 

One final note, a s depends very much on details of 
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FIG. 24: (Color online) Surface coefficient of the level-density 
parameter obtained for the A ~170 systems as a function of 
excitation energy. The dotted line indicates the value ob- 
tained from the fit in Fig. 1171 Other values of the surface 
coefficient, from theprescriptions of Toke and Swiatecki[44l| 
and Ignatyuk et aZ.|Ig], are also indicated. 



how the mean-field potential changes with deformation. 
IIV Al the diffuseness parameter d of the nuclear 



In Sec. 

potential was made angle dependent (Eq. I26fl so that the 
diffuseness perpendicular to the nuclear surface is con- 
stant. The parameter d is actually the diffuseness along 
the radial direction. If instead d is set to be constant 
as in many studies, then the extracted a s values were 
found to be negative! In this case, the mean diffuseness 
perpendicular to the nuclear surface decreases with defor- 
mation. As the level-density parameter is quite sensitive 
to the diffuseness, this leads to the calculated reduction 
of the level-density parameter with deformation. Thus 
for large deformations, it is important to make d angle 
dependent. 
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VIII. DISCUSSION 

In the preceding section, it was shown that the 
two methods for calculating the contribution from the 
positive-energy single-particles levels give similar results, 
apart from the deformation dependence. The difference 
between the two methods are the inclusion of both the 
negative background and the wide resonances in the sub- 
traction method. These two contributions have opposite 
effects on the level density, and thus it seems they par- 
tially cancel each other. Also for heavy systems with 
large Coulomb barriers, the low, positive-energy proton 
states are all narrow resonances and thus give almost 
identical results for the two methods. Thus the proton- 
rich side of the chart of nuclides is less sensitive than 
the neutron-rich side to the continuum, at least for the 
heavier systems. 

If the Gamov method is considered preferable, then 
consideration must be given to the value of T$. Possibly 
To is related to the CN lifetime, i.e., only Gamov levels 
of lifetime greater than the CN should be considered. As 
the CN lifetime decreases with excitation energy, then Tq 
should increase. An increasing value of To with excitation 
energy would lead to a reduction in the predicted strong 
peaking of the level density at sphericity. Another possi- 
bility is that Tq should be related to the average spread- 
ing width of levels near the Fermi surface. This would 
lead to a much smaller excitation-energy dependence of 

r . 

It is important to remember that the results of this 
work relate to how the treatment of positive-energy 
single-particle levels affects the level density in the 
independent-particle model. Apart from pairing, the cal- 
culations did not include any other many-body effects. 
Also they did not allow for selfconsistancy between the 
assumed nuclear potential and predicted density distri- 
butions of the nucleons. These effects may to lead to 
deviations from the predicted behavior. In fact, exper- 
imentally it is known that the level density for hot Yb 
nuclei cannot be described by a backshifted Fermi-gas ex- 
pression, but an important excitation-energy dependence 
of the level-density parameter is neededpj. The latter is 
consistent with the variation of the frequency-dependent 
effective nucleon mass with temperature|50j . 

Some of the effects ignored in this work may influence 
the n-p asymmetry dependence. For example, if the dif- 
fuseness parameter of the nuclear potential increases for 
nuclei close to the drip lines, then this will enhance the 
level-density parameter for these systems. Also differ- 
ences between the neutron and proton effective masses 
and their dependence on asymmetry may also be impor- 
tant. 

Al-Quraishi et al. 3 had suggested that a restriction of 
the positive-energy states to narrow resonances, as in the 
Gamov method, would lead to an important n-p asym- 
metry dependence of the level-density parameter. They 
fit the density of known levels for 20< A <70 with the 



form 



a c = ai^lexp |a 2 [Z - Zp (A)] 2 j 



(65) 



where Zp (A) is the proton number of the /3-stable nu- 
cleus of nucleon number A and a\ and are the fit 
parameters. This is called case C in Ref. 0. The dashed 
curves in Figs. 116b and 116b show ac for A~60 and 40, 
respectively. If this strong n-p asymmetry dependence is 
real, the calculations of this work suggest it cannot be ex- 
plained by the treatment of the continuum as assumed in 
the justification of Eq.|HHl Al-Quraishi et al. also consid- 
ered another fit (case B) based on isospin considerations: 



a-B = 03^4 exp 04 (N — Z) 



(66) 



where now 03 and 04 are the fit parameters. The dot- 
dashed curves in Figs. IHSb and iTob show the resulting 
values of as- In this case, the n-p asymmetry depen- 
dence is not as strong as in case C, but still stronger 
than our calculations for A=60. The large difference be- 
tween cases B and C suggests the fits do not constrain the 
level-density parameter for very neutron and proton-rich 
nuclei. In Ref. 

El. 

no significant asymmetry dependence 
of the level-density parameter was observed for 152 Yb and 
160 Yb CN with excitation energy greater than 100 MeV. 
This is consistent with the dependence calculated with 
both methods of this work (see Figs. 116b). 

It is also important to note again that the small n- 
p asymmetry dependence observed in these calculations 
is only true for T < E™™ t . Extending the calculations 
above T = E™ 1 ^ leads to a much larger dependence, al- 
though this region is not relevant to the statistical model. 
In Ref. [18j is was noted that for calculations where the 
number of single-particle levels is finite, such as in the 
Gamov method, then there is a maximum excitation en- 
ergy of the CN. The level density as a function of exci- 
tation energy must peak and then approach zero at this 
maximum excitation energy. The peak value of the level 
density corresponds to T=oo and higher excitation ener- 
gies have negative temperatures. The details of this be- 
havior would be very dependent on the n-p asymmetry. 
Although a negative temperature may be appropriate for 
the CN, it is certainly not meaningful for the gas. Thus 
our model of CN in equilibrium with the surrounding gas 
breaks down. However, it is not clear that this is of any 
relevance for the statistical model. 



IX. CONCLUSIONS 

The effects of continuum positive-energy neutron and 
proton levels on the nuclear level density has been in- 
vestigated. The use of the independent-particle model 
allowed for a broad survey of how these continuum cor- 
rections modify the level density over the entire chart of 
nuclides. Two methods for calculating the contributions 
of these positive-energy levels were investigated. In the 
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subtraction method, the single-particle level density is 
determined from the scattering phase shifts. The result- 
ing single-particle level density has contributions from 
narrow and wide resonances and a negative background. 
In the Gamov method, the single-particle level density 
is calculated from the Gamov states and only the con- 
tributions from the narrow resonances are considered. 
From the bound states and these two prescriptions for 
the positive-energy states, the entropy and level density 
are calculated as a function of temperature and excitation 
energy. These calculations ignored all many-body effects 
apart from the pairing interaction. At large excitation 
energies where shell effects melt, the level density fol- 
lowed a backshifted Fermi-gas expression. Also, the de- 
duced level-density parameters were quite similar for the 
two methods. They depend on A with very little depen- 
dence on the n-p asymmetry of the nucleus. The biggest 
asymmetry dependence was for the very exotic systems 
near the drip lines where a small reduction in the level- 
density parameter was sometimes found. For the heavier 
systems, the prescription of Ignatyuk et q/.|45j| accounted 
for the variation in level density at low excitation energies 
where shell effects are still important. 



The largest differences arising from the use of the two 
methods was the predicted deformation dependence of 
the level density. At high excitation energies, the Gamov 
method predicted the level density peaked strongly for 
spherical systems whereas in the subtraction method the 
deformation dependence was rather flat near sphericity. 
This suppression in the relative level densities of de- 
formed to spherical systems in the Gamov method would 
lead to a reduction in the predicted fission width and 
may help explain the large prescission light-particle mul- 
tiplicities observed in fission reactions |48|. 
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